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Aromatic  cations  such  as  imidazolium  and  its  derivatives  have 
delocalized  charges  and,  therefore,  are  often  used  to  synthesize 
functional  nanomaterials,  low-melting  salts,  and  ionic  liquids.1 
Recently,  energetic  low-melting  salts  and  ionic  liquids  based  on 
triazolium  and  its  derivatives  have  received  considerable  attention.2’3 
A  characteristic  property  of  aromatic  compounds  is  their  ability  to 
form  7t— n  stacking  structures,  which  can  have  significant  influence 
on  properties.  Cation— cation  n— n  stacking  structures  are  less 
common  but  have  been  observed  in  organometallic  crystals4  and 
imidazolium  salts.5  Cation— cation  n—n  stacking  dimers  can  be 
formed  in  a  dielectric  solvent  due  to  significant  screening  of  the 
charge— charge  repulsion  forces.6  However,  such  species  in  small 
clusters  have  not  been  reported.  Calculations  reported  in  this  work 
suggest  that  cation— cation  n—n  stacking  structures  can  be  formed 
between  two  1,2,4-triazolium  cations  in  small  clusters  of  the  two 
cations  and  two  anions. 

The  ab  initio  calculations  were  performed  with  the  electronic 
structure  package  GAMESS.7  Geometries  were  optimized  using 
second-order  perturbation  theory  (MP2)8  with  the  aug-cc-pVDZ9 
basis  set.  Single  point  energies  were  calculated  for  selected 
structures  using  MP2  with  the  cc-pVDZ,  aug-cc-pVDZ,  aug-cc- 
pVTZ,  and  aug-cc-pVQZ  basis  sets9  and  coupled  cluster  singles, 
doubles  with  perturbative  triples  [CCSD(T)]10  with  the  cc-pVDZ, 
and  aug-cc-pVDZ  basis  sets. 

Numerous  structures  can  be  found  for  two  pairs  of  1,2,4- 
triazolium  and  dinitramide  ions  or  two  pairs  of  1,2,4-triazole  and 
dinitramine  neutral  molecules.  An  extensive  geometry  search  was 
performed  for  the  ionic  and  neutral  structures.  The  lowest  energy 
MP2/aug-cc-pVDZ  structures  are  shown  in  Figure  1.  In  the  most 
stable  ionic  two-pair  structure,  each  1,2,4-triazolium  forms  two 
hydrogen  bonds,  via  the  H’s  on  the  N  atoms,  to  the  O  atoms  of  the 
dinitramide  ions.  Interestingly,  this  structure  exhibits  Jt—n  stacking 
of  the  two  cationic  1,2,4-triazolium  rings.  The  interplane  distance 
is  ~3.2  A,  with  a  parallel  displacement  of  ~  1.4  A.  The  nitro  groups 
of  the  dinitramide  ions  are  approximately  perpendicular  to  each 
other.  In  the  most  stable  neutral  two-pair  structure,  each  1,2,4- 
triazole  forms  one  N— H— N  hydrogen  bond  to  a  dinitramine.  There 
is  Ji—Jt  stacking  between  the  two  approximately  parallel  1,2,4- 
triazole  rings,  with  an  interplane  distance  of  ~3.2  A  and  a  parallel 
displacement  of  ~1.4  A.  For  comparison,  an  additional  ionic 
structure  was  created  from  the  neutral  structure  by  repositioning 
protons  onto  the  1,2,4-triazole  molecules.  Upon  geometry  optimiza¬ 
tion,  the  protons  transfer  back  to  the  dinitramide  ions  to  form  neutral 
pairs  with  no  energy  barrier. 

It  is  of  interest  to  determine  the  cluster  size  at  which  the  ion 
pair  structures  become  more  stable  than  the  corresponding  neutral 
pair  structures.  A  previous  study  predicted  that  ion  pair  dimers  are 
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Figure  1.  MP2/aug-cc-pVDZ  optimized  structures  of  two  pairs  of  1,2,4- 
triazolium  (1,2,4-triazole)  and  dinitramide  (dinitramine)  molecules.  H  is 
white,  C  is  gray,  O  is  red,  N  is  blue. 


typically  higher  in  energy  than  neutral  pair  dimers.3  The  MP2/aug- 
cc-pVDZ  energy  of  the  ionic  tetramer  in  Figure  1  is  1.5  kcal/mol 
lower  than  that  of  the  neutral  one.  Including  zero  point  energies 
(ZPE)  obtained  with  MP2/aug-cc-pVDZ  harmonic  vibrational 
frequencies,  the  ionic  tetramer  is  1.2  kcal/mol  lower  in  energy  than 
that  of  the  neutral  one.  The  MP2  method  tends  to  predict  higher 
energies  for  ionic  species  versus  neutral  species.3  Since  the 
computational  cost  of  CCSD(T)/aug-cc-pVDZ  energy  calculations 
of  these  two  tetramer  structures  is  prohibitive,  these  energies  were 
approximated  from  the  MP2/aug-cc-pVDZ  energies  by  estimating 
the  electron  correlation  energy  differences  using  three  independent 
methods:  (1)  the  differences  between  the  MP2/cc-pVDZ  and 
CCSD(T)/cc-pVDZ  energies  of  the  tetramers,  (2)  the  differences 
between  the  MP2/aug-cc-pVDZ  and  CCSD(T)/aug-cc-pVDZ  ener¬ 
gies  of  the  12  pairs  of  dimers  in  these  two  tetramers,  and  (3)  the 
differences  between  the  MP2/aug-cc-pVDZ  and  CCSD(T)/aug-cc- 
pVDZ  energies  of  the  eight  monomers  in  these  two  tetramers.  Using 
these  three  methods,  the  estimated  CCSD(T)/aug-cc-pVDZ  energy 
of  the  ionic  tetramer  is  lower  than  that  of  the  neutral  tetramer  by 
6.0,  7.6,  and  8.0  kcal/mol,  respectively.  Including  ZPE,  these  three 
values  become  5.7,  7.3,  and  7.7  kcal/mol,  respectively.  Further 
corrections  for  the  basis  set  effects  from  aug-cc-pVDZ  to  aug-cc- 
pVTZ  and  to  aug-cc-pVQZ  are  small.  One  can  therefore  conclude 
that,  for  two  pairs  of  1,2,4-triazolium  and  dinitramide,  ionic 
structures  are  more  stable  than  the  corresponding  neutral  structures. 

The  driving  force  for  the  formation  of  the  unusual  cation— cation 
n— n  stacking  structure  is  also  very  interesting.  Geometry  optimiza¬ 
tion  of  the  ionic  structure  using  HF,  which  does  not  contain  electron 
correlation,  results  in  an  open-book  structure  in  which  the  two  1,2,4- 
triazolium  planes  form  a  dihedral  angle  of  ~120°  (no  Jt—Ji 
stacking).  Clearly,  electron  correlation  energy  (a  major  part  of  the 
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Figure  2.  MP2/aug-cc-pVDZ  optimized  structures  of  two  pairs  of  1,2,4- 
triazolium  chloride  ions.  H  is  white,  C  is  gray,  N  is  blue.  Cl  is  green. 

aromatic  ji—ji  attraction)  in  the  MP2  method  is  the  driving  force 
for  the  formation  of  the  cation— cation  jt—jt  stacking  structure  from 
the  open-book  non-stacking  structure. 

The  tetramer  interaction  energy  and  the  six  isolated  dimer 
interaction  energies  of  the  four  monomers  in  the  ionic  tetramer 
structure  in  Figure  1  were  derived  by  taking  the  corresponding 
energy  differences  between  the  tetramer  and  isolated  dimers  and 
monomers,  with  no  correction  for  basis  set  superposition  errors. 
The  same  geometries  as  in  the  tetramer  were  used  for  the  isolated 
dimers  and  monomers,  and  the  aug-cc-pVDZ  basis  set  was  used. 
Using  HF,  the  tetramer  interaction  energy  is  —217.1  kcal/mol;  the 
repulsion  energies  between  the  two  1,2,4-triazolium  cations  and 
the  two  dinitramide  anions  are  +77.1  and  +52.9  kcal/mol, 
respectively;  the  attraction  energies  between  the  two  cations  and 
two  anions  are  —96.8,  —89.3,  —89.3,  and  —96.8  kcal/mol  (sym¬ 
metric  due  to  the  symmetry  of  the  tetramer  structure).  The  —217.1 
kcal/mol  tetramer  interaction  energy  is  roughly  the  sum  (—242.2 
kcal/mol)  of  the  dimer  interactions,  with  a  +25.1  kcal/mol 
difference  due  to  many-body  effects.  These  HF  interaction  energies 
are  dominated  by  Coulomb  interactions.  For  example,  a  reduced 
variational  space* 1 2 3 4 5 6 7 8 9 10 11  analysis  of  the  two  1,2,4-triazolium  cations  leads 
to  the  following  terms  (in  kcal/mol):  total  =  +77.1,  Coulomb  = 
+76.9,  exchange-repulsion  =  +4.9,  polarization  =  —3.8,  charge 
transfer  =  —0.5,  and  mixing  terms  =  —0.4.  In  general,  the  driving 
force  for  the  formation  of  ionic  clusters  is  the  Coulomb  attraction 
between  the  opposite-charge  ions,  which  is  typically  much  stronger 
than  the  repulsion  between  the  same-charge  ions. 

However,  as  shown  above,  additional  dispersion  attraction  due 
to  electron  correlation  is  necessary  for  the  formation  of  the  cationic 
ji—ji  stacking  structure  from  the  (HF)  open-book  non-stacking 
structure.  The  MP2  tetramer  interaction  energy  is  —266.7  kcal/ 
mol,  49.6  kcal/mol  stronger  than  the  HF  one,  due  to  the  electron 
correlation  attraction;  the  MP2  interaction  is  roughly  the  sum 
(—294.5  kcal/mol)  of  the  MP2  dimer  interactions,  differing  by 
+27.8  kcal/mol  due  to  many-body  effects.  The  CCSD(T)  tetramer 
interaction  energy  is  not  available;  the  sum  of  the  CCSD(T)  dimer 
interactions  is  —290.9  kcal/mol,  very  close  to  the  MP2  one.  For 
the  two  1,2,4-triazolium  cations,  the  HF,  MP2,  and  CCSD(T) 
interaction  energies  are  +77.1,  +69.3,  and  +70.4  kcal/mol, 
respectively;  the  electron  correlation  corrections  to  the  HF  interac¬ 
tion  energy  predicted  by  MP2  and  CCSD(T)  are  —7.8  and  —6.7 
kcal/mol,  respectively. 

On  the  basis  of  the  structure  of  two  1,2,4-triazolium— dinitramide 
ion  pairs  (Figure  1),  an  initial  structure  of  two  1,2,4-triazolium— 
chloride  ion  pairs  was  constructed  and  then  optimized  using  MP2/ 
aug-cc-pVDZ.  In  the  optimized  structure  shown  in  Figure  2,  each 
1,2,4-triazolium  ion  forms  two  N— H— Cl  hydrogen  bonds.  Appar¬ 


ently,  ji—ji  stacking  of  the  two  cationic  1,2,4-triazolium  rings  is 
formed  in  this  structure  in  spite  of  the  small  sizes  of  the  chloride 
ions.  The  interplane  distance  is  ~2.9  A,  with  a  parallel  displacement 
of  ~  1.0  A.  The  structure  is  not  symmetric  due  to  the  parallel  up/ 
down  displacement  (front  view  in  Figure  2).  Clearly,  very  strong 
N— H— Cl  hydrogen  bonds  are  formed  in  this  structure  so  that  the 
four  H  atoms  bend  inward  (The  N---C1  distances  are  2.910  and 
2.929  A).  A  neutral  structure  of  two  pairs  of  1,2,4-triazole  and 
hydrogen  chloride  (HC1)  molecules  was  constructed  starting  from 
the  1,2,4-triazole— dinitramine  neutral  tetramer  in  Figure  1  and  then 
optimized  with  MP2/aug-cc-pVDZ.  The  neutral  structure  optimized 
to  the  ionic  structure  shown  in  Figure  2,  with  no  energy  barrier. 
Again,  electron  correlation  attraction  is  the  driving  force  for  the 
formation  of  the  7t— tt  stacking  structure:  HF  geometry  optimization 
of  the  structure  leads  to  a  planar  structure  in  which  all  the  atoms 
are  exactly  in  the  same  plane. 

In  conclusion,  ab  initio  calculations  suggest  that  cation— cation 
tt— 71  stacking  structures  can  exist  in  very  small  ionic  clusters  such 
as  two  1,2,4-triazolium  cations  and  two  dinitramide  or  chloride 
anions.  The  structure  motifs  and  interaction  patterns  provide  new 
understanding  of  ionic  materials  with  aromatic  rings. 

Acknowledgment.  This  work  was  supported  by  a  grant  from 
the  Air  Force  Office  of  Scientific  Research.  The  calculations  were 
performed  at  the  MSCF  high-performance  supercomputer  center 
at  Pacific  Northwest  National  Laboratory,  at  DoD  Aeronautical 
Systems  Center,  Wright-Patterson  Air  Force  Base,  OH,  and  the 
Naval  Oceanographic  Office,  John  Stennis  Space  Center,  MS,  under 
a  DoD  Challenge  grant,  and  on  a  computer  obtained  with  startup 
funds  from  the  University  of  Nebraska— Lincoln. 


References 

(1)  Wilkes,  J.  S.;  Zaworotko,  M.  J.  J.  Chem.  Soc.,  Chem.  Commun.  1992, 
965-967. 

(2)  (a)  Drake,  G.;  Hawkins,  T.;  Brand,  A.;  Hall,  L.;  McKay,  M.;  Vij,  A.; 
Ismail,  I.  Propellants  Explosives  Pyrotechnics  2003,  28,  174—180.  (b) 
Singh,  R.  P.;  Verma,  R.  D.;  Meshri,  D.  T.;  Shreeve,  J.  M.  Angew.  Chem., 
Int.  Ed.  2006,  45,  3584-3601. 

(3)  Schmidt,  M.  W.;  Gordon,  M.  S.;  Boatz,  J.  A.  J.  Phys.  Chem.  A  2005. 
109,  7285-7295. 

(4)  Janiak,  C.  J.  Chem.  Soc.,  Dalton  Trans.  2000,  3885-3896. 

(5)  (a)  Wilkes,  J.  S.;  Zaworotko,  M.  J.  Supramol.  Chem.  1993,  191-193.  (b) 
Bujak,  M.;  Zaleski,  J.  Z.  Naturforsch.  2002,  57,  157—164.  (c)  Sieron,  L. 
Acta  Crystallogr.  Sect.  E  2007,  63,  0904-0905. 

(6)  Buisine,  E.;  de  Villiers,  K.;  Egan,  T.  J.;  Biot,  C.  J.  Am.  Chem.  Soc.  2006, 
128,  12122-12128. 

(7)  (a)  Schmidt,  M.  W.;  Baldridge,  K.  K.;  Boatz,  J.  A.;  Elbert,  S.  T.;  Gordon, 
M.  S.;  Jensen,  J.  H.;  Koseki,  S.;  Matsunaga,  N.;  Nguyen,  K.  A.;  Su,  S.  J.; 
Windus,  T.  L.;  Dupuis,  M.;  Montgomery,  J.  A.  J.  Comput.  Chem.  1993, 
14,  1347—1363.  (b)  Gordon,  M.  S.;  Schmidt,  M.  W.  Advances  in  electronic 
structure  theory:  GAMES S  a  decade  later.  In  Theory  and  applications  of 
computational  chemistry,  Dykstra,  C.  E.,  Frenking,  G.,  Kim,  K.  S., 
Scuseria,  G.  E.,  Eds.;  Elsevier:  Amsterdam,  2005. 

(8)  (a)  M0ller,  C.;  Plesset,  M.  S.  Phys.  Rev.  1934,  46,  618-622.  (b)  Pople, 

J.  A.;  Binkley,  J.  S.;  Seeger,  R.  Int.  J.  Quantum  Chem.  1976,  1  —  19.  (c) 
Aikens,  C.  M.;  Webb,  S.  P.;  Bell,  R.  L.;  Fletcher,  G.  D.;  Schmidt,  M. 
W.;  Gordon,  M.  S.  Theor.  Chem.  Acc.  2003,  110,  233-253. 

(9)  Dunning,  T.  H.  J.  Chem.  Phys.  1989,  90,  1007-1023. 

(10)  (a)  Raghavachari,  K.;  Trucks,  G.  W.;  Pople,  J.  A.;  Head-Gordon,  M.  Chem. 
Phys.  Lett.  1989,  479—483.  (b)  Piecuch,  P.;  Kucharski,  S.  A.;  Kowalski, 

K. ;  Musial,  M.  Comput.  Phys.  Commun.  2002,  149,  71—96.  (c)  Bentz,  J. 

L. ;  Olson,  R.  M.;  Gordon,  M.  S.;  Schmidt,  M.  W.;  Kendall,  R.  A.  Comput. 
Phys.  Commun.  2007,  176,  589—600.  (d)  Olson,  R.  M.;  Bentz,  J.  L.; 
Kendall,  R.  A.;  Schmidt,  M.  W.;  Gordon,  M.  S.  J.  Chem.  Theory  Comput. 
2007,  3,  1312-1328. 

(11)  Stevens,  W.  J.;  Fink,  W.  H.  Chem.  Phys.  Lett.  1987,  139,  15-22. 
JA076406U 


J.  AM.  CHEM.  SOC.  ■  VOL.  130,  NO.  2,  2008  393 


